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Abstract 



oo This paper presents a detailed numerical study of the effect of 
On 

C\ focusing on the conversion efficiency of low-loss singly- 

Ojt) resonant parametric oscillators with collinear focusing of 

^ pump and signal. Results are given for the maximal pump 

G\ depletion and for pump power levels required for various 

amounts of depletion, as functions of pump and signal confocal 

^ parameters, for ]q/kp=0.33 and 0.50. It is found that the ratio 

O 

of pump depletion to maximal depletion as a function of the 
ratio of pump power to threshold power agrees with the plane- 
wave prediction to within 5%, for a wide range of focusing 

o 

conditions. The observed trends are explained as resulting 
j£ from intensity and phase dependent mechanisms. 



Introduction 

The recent demonstrations of continuous-wave operation for singly-resonant 
optical parametric oscillators (SRO) with significant pump depletion in collinear 
focusing configurations [1-6] have stimulated renewed interest in these devices, and call 
for a more detailed analysis of diffractive effects. An exact calculation of thresholds for 
collinear focusing was presented in [7]. Theoretical studies of SROs above threshold so 
far have used the plane-wave [2,4,8] or lowest-order Gaussian mode [1,2,4,5,9] 
approximations. A numerical model of a pulsed SRO with strong birefringent walk-off 
is found in reference [10]. In [11], the author studied the effects of residual dispersion 



and diffraction on threshold, pump depletion and transverse mode content in a model of 
a low-loss SRO with collinear focusing at the center of the crystal, and a confocal 
parameter of signal and pump equal to the crystal length. It was found that the 
transverse mode content of the idler depends strongly on the residual dispersion, but 
only weakly on the amount of pump depletion. At optimal dispersion, about 77% of the 
total idler power is in the zero order radial mode, 22% in the first mode, and 1% in 
higher modes. Maximal pump depletion of 95.5% occurs at a pump power 2.52 times 
above threshold. 

In this paper, we study the effect of various degrees of focusing for pump and 
signal on threshold and pump depletion for the same model. In section 1, we present the 
coupled wave-equations, define intrinsic units for the length scales and power, and 
discuss the numerical solution procedure. In section 2, we present our results for the 
maximal degree of pump depletion and for the required pump power for threshold and 
various degrees of pump depletion as functions of the confocal parameters of pump and 
signal. In section 3, we discuss the physical mechanisms responsible for the observed 
numerical trends. In section 4, we discuss the limitations of our model, and summarize 
our findings. 

1. Field equations, resonator modeling and numerical procedures 
a. Field equations, intrinsic units and photon conservation 

The electric field Ef and the total power Pf of the pump, signal and idler fields 
(F=P,S,I) are given in terms of the complex field amplitude Af as: 

EF=real(Apexp(i<£F)X <X>F=kFZ-iC0t, kF=COFnF/c=27tnF/^F (1) 




For steady-state continuous-wave operation, the coupled wave equations for the 
amplitudes Af in a second-order non-linear medium are: 



A F?Z = (i/2k F )(A F ,xx + A F , yy ) + Qf 



(2) 



Qp=i Kpexp(-iAz)AiAs, Qs=i Ks exp(iAz)ApAi*, Qi=i Ki exp(iAz)ApAs* , 
KF=C0Fd/nFC=27id/nF^F, A=kp-ks-ki-2jt/l. 

Here d is the non-linear coefficient, and A=kp-ks-ki-27t/l is the residual 

dispersion for a periodically poled crystal with poling period 1. The subscripts x,y and z 
denote partial differentiation with respect to the corresponding spatial coordinates. For 
example, for periodically poled LiNb03 (PPLN) at ^p=1.06 |im, one has approximately 

d=14 and n=2.2. To facilitate the application of our numerical data to various 
experimental situations, we introduce arbitrary units of propagation distance zo and 

wave-vector ko, which define corresponding units ro, Apo and Pfo for the radial 

distance, complex amplitude and power as follows: 

ro 2 =z /ko , A P o=1/(zoVk S Kt), Aso=1/(z(Wk P Ki), Aio=1/(z Vk S k P ) (3) 

With these units, the wave equations 2 and the expression for the beam power Pf 
take the form: 




(i/2k'p)(A' P , x 



x' + A'p, 



y'yO + i exp(-iA'z')A'iA's , 



(4.a) 




(i/2k'i)( A'i,x'x' + A' 



I, y y) + i exp( iA'z')A' P A's* , 



(4.b) 




(i/2k's)(A' S ,x'x' + 



A'S, 



y-yO + i exp( iA'z')A' P A'i* , 



(4.c) 



z=z'zq, kF=kF'ko, x=x'ro, y=y'ro, AF=aF'AF0, A=A'/zo . 




(5) 



P0P=£0 c np ns ni Xs ^1/(871 d 2 zq ko). 



The units Pqf of power are related by: 



X? Pop = Xi Poi = Pos- 



(6) 



The dimensionless quantities rp are proportional to the total rates of photons 

passing through the crystal in unit time. They satisfy the following relationships, which 
are expressions of photon number conservation: 



In our following analysis, we set zo=L/2, where L is the crystal length, and 
ko=kp. In this case, the unit Pop of pump power is given by: 

P0P=£0 c n S ni Xs h ^p/(87t 2 d 2 L) = 3.352 10 4 n S ni Xs h Xp/(d 2 L) , (8) 
[X] = jam , [d] = pm/V , [L] = mm, [Pop] = Watt . 

The expressions in square brackets denote the units for which the above 
numerical expression is valid. The unit Pop of pump power is related to the quantity K3 

of reference [7] by: 



b. Resonator modeling 

We now use the photon conservation equations 7 to derive relationships between 
resonator input and output in a SRO with a low-loss ring resonator. We denote the 
photon rates at the input and output faces of the crystal by rFin and rF ut- Then the 
following relationships hold: 



d(rp +ri)/dz' = d(rp + rs)/dz' = d(rs - ri)/dz' = 0. 



(7) 



Pop = 8/((l + X s /Xj>) K 3 ). 



(9) 



riin =0, rsin = rs, 



rpin = rp, 



(10) 



riout = n, rSout = rs + rr, 



rpout = rp . n . 



If T«l is the total resonator loss for the signal, one has 



rSin = (1-T) rsout = rs + n - Trs, n = Trs « rs. 



(11) 



In this case, the total signal power in the cavity stays almost constant. If the input 
pump beam is a zero order radial mode (Gaussian), the zero order radial signal mode 
has lowest threshold. We found that the threshold for the first radial signal mode is 
about four times the threshold for the zero order mode. Therefore, it is a valid 
approximation to replace the signal amplitude in equations 4.a and 4.b by the amplitude 
corresponding to a zero order mode with constant total power. Equations 4. a and 4.b 
are then linear in the pump and idler amplitude, and the total idler output power is a 
linear function of pump input power, with a conversion coefficient D(rs) depending on 

the signal power: 



From equations 1 1 and 12, the pump depletion Dp and pump power are obtained 
as functions of rs as: 



The normalized pump power Rp as a function of rs depends only on the focusing 

geometry and ratio of idler to pump wave length, while all dependence on the resonator 
loss and the non linear coefficient is contained in T and Pop. In the limit of rs -> 0, 

equation 13.b gives the normalized threshold pump power Rthr-The gain coefficient 

D(rs) is calculated numerically by the procedure described below. 

c. Numerical algorithms and parameters 

Equations 4.a and 4.b were integrated numerically by a split step algorithm, with 



ri = D(rs) rp. 



(12) 



Dp = n/rp = D(rs), 



(13.a) 



rp= ri/D(rs) = T rs/D(rs) = T Rp , Pp=rpPop=RpTPop. 



(13.b) 



an explicit half- step of parametric conversion followed by a full step of diffraction and 
an implicit half step of conversion. The diffraction step uses the Cayley transform of a 
discrete version of the radial Laplacian which conserves photon number. This 
procedure is unconditionally stable [12]. The use of the radial Laplacian restricts us to 
radially symmetric beams, which is a reasonable approximation for many experimental 
situations. It results in a considerable reduction in computation time as compared to a 
full two-dimensional transverse integration. 

As discussed above, the signal beam was set to a Gaussian with constant power rs 

throughout the crystal. The amplitude of the pump beam at the crystal entry face was 
set to be that of a Gaussian with unit power rp = 1.0, while the amplitude of the idler 

beam was set to zero. The total output power of the idler at the crystal exit face was 
calculated for a sequence of values rs of signal power, to give the conversion or pump 

depletion coefficient D(rs), according to equation 12. At each value of rs, the residual 

dispersion parameter A was adjusted to give maximal idler output power. Equations 
13. a and 13.b were then used to obtain the pump depletion D(Rp) as a function of the 

actual normalized pump input power Rp corresponding to the specified values of signal 

power rs, according to the resonator condition, equation 11. 

We verified our numerical procedure extrinsic by propagating zero and first 
order radial modes without conversion, and plane wave beams with full conversion. 
Agreement of the numerical beam amplitudes with the known analytical solutions was 
better than 0. 1 % of the maximal amplitudes, for the choice of radial and longitudinal 
step sizes given below. Larger errors occurred for confocal parameters smaller than 
about one half the crystal length at radial distances where the beam amplitudes had 
decreased to less than about 5% of their maximal values and longitudinal distances more 
than a Raleigh range from focus. These parts of the beam profile contribute less than 
0.25% of the total beam power, and the errors were therefore tolerated. In [11], we 
also found very close agreement between the numerical results obtained from the 
procedure described here, and a numerical integration of coupled mode equations. 
Furthermore, the threshold power calculated was found to agree with the results of 
reference [7]. 

Intrinsic verification was obtained by reducing the radial and longitudinal step 
sizes by one half, resulting in changes of less than 0.2% of the calculated idler power, 
which is consistent with the results from the extrinsic verification described above. The 



crystal length L was defined to be two units of length: L=2Lo , and was divide into 80 
integration steps. For the smallest pump confocal parameter bp=0.3L considered here, 
this gives 24 steps for the region within on Raleigh range from the focus. With a unit 
ko=kp of wave vector, a pump beam with bp=L=2Lo has a waist radius of w=1.4wo 

units wo of radial distance. We chose the maximal radial distance to be 14 units, and the 

number of radial integration steps to be 200. This gives 20 steps within the waist radius 
for bp=L, and more than 10 steps for bp=0.3L. 

The signal power rs was increased from rs =0.0001 by 0.2 to a final value of 

4.0001. Maximal conversion was found in the range of rs=2.8 to 3.2. From the 

resulting discrete series of pump depletion or conversion coefficients the maximal 
conversion and pump power levels at fixed conversion were determined by 
interpolation. 

2. Numerical results for pump depletion and pump power levels 

Figures 1 to 9 summarize results for various parameters, which characterize the 
conversion efficiency of the SRO as functions of the inverse focusing parameter 
zp'=zp/Lo=bp/L. Figures 1 to 6 present results for fixed signal confocal parameter 

zs -bs/L, with figures 1 to 3 for the case lq/kp=0.33 and figures 4 to 6 for ki/kp=0.50. 

Figures 7 to 9 give corresponding results for the case where the pump and signal 
confocal parameters are equal: zp'=zs'. 

a. Results for ki/kp=0.33, bs/L=1.0, 0.7 and 0.5. 

Figure 1 shows the maximal pump depletion D max as a function of zp', for 

zs -1.0, 0.7 and 0.5. The largest value is always found at zp' smaller than or equal to 

zs'. For example, at zs -1.0, D max peaks between zp'=0.8 and 0.9 at about 0.965, while 

for zs -0.5, the peak value of 0.977 is found at zp'=0.5. Thus, for maximal conversion, 

it is necessary to focus the pump not less tightly than the signal. The penalty for 
violating this requirement is not severe, however. For zs -1.0, one stays within 5% of 

the peak value, that is above 0.91, for zp' in the range from 0.5 to 1.5. 

Figure 2 shows the pump power levels Rp at which maximal depletion D max and 



95% of maximal depletion, D=0.95D max , are achieved. For zs -1.0, the lowest pump 

level required is Rp=3.12 for maximal depletion, and Rp=2.34 for 95% of maximal 

depletion. The corresponding value of zp' is between 0.7 and 0.8, somewhat below the 

value at which D max peaks. For zs -0.5, the minimal levels are lower: Rp=2.87 and 

2.12, and they occur for the same zp'=0.5 as for peak D max . 

Figures 3. a to 3.c show the pump levels Rp required to reach threshold D=0.0, 

and the values D=0.2, 0.4, 0.6 and 0.8. The threshold values agree with those published 
in [7] (see also [10]). Minimal power levels are found in the same range for zp' as in 

figure 2: zp'=0.7 to 0.8 for zs -1.0 (figure 3. a), zp'=0.6 for zs -0.7 (figure 3.b) and 
zp'=0.5 for zs -0.5 (figure 3.c). 

Figure 3.d finally gives the ratio Pp/Pthr=Rp/Rthr of pump power at D=0.95D max 
(figure 2) to threshold power (figures 3. a to 3.c). It is very interesting to note that this 
ratio falls into a narrow range from 1.99 to 2.07, a variation of only 4%. Thus, as a 
rule of thumb for a wide range of focusing conditions one expects D=0.95D max to 

occur at about 2 times above threshold. Violation of this rule will indicate the presence 
of nonlinear effects other than parametric conversion. 

b. Results for ki/k P =0.50, bs/L=1.0, 0.7 and 0.5. 

Figures 4 to 6 present corresponding results for the degenerate case 
ki/kp=ks/kp=0.50. The maximal depletion D max is somewhat higher than for 

kp/kp=0.33, although not by a large amount. For zs -0.5, it peaks just below D max =0.98 

at zp'=0.5. 

The pump levels are lower than for kp/kp=0.33. For example, at zs -0.5, the 

minimal level for D=0.95D max is Rp=2.18, while it is Rp=2.34 for ki/kp=0.33. The 

same trend is found for the relative pump power levels P/Pthr at D=0.95D max , shown in 

figure 6.d. These lie in the range from 1.93 to 1.99, as compared to the range 1.99 to 
2.07 for ki/k P =0.33. 



c. Results for bs=bp, ki/kp=0.33, 0.50. 



Figures 7 to 9 present results for the case where signal and pump have equal 
confocal parameter. Peak D max and lowest power levels are found for tighter focusing, 

around zp'=zs -0.4. The degenerate case ki/kp=0.50 requires significantly lower pump 

power levels for given depletion D than ki/kp=0.33 at tighter focusing zp'=zs'<1.0, 

while for larger confocal parameters the difference is less. Differences in Dmax 
between the two cases are less than 1% for all values of zp'=zs'. 

d. Relative pump depletion as a function of relative pump power. 

Figure 10 gives the relative depletion D/D max as a function of relative pump 

power P/Pthr, showing one near optimal case each for ki/kp=0.33 and 0.50, and the 

result of the plane wave theory. The curve for kp/kp=0.50 follows the plane wave curve 

more closely, staying below the plane wave curve until Dmax is reached, and declining 
more gradually afterwards. This behavior is more pronounced for ki/kp=0.33. 

All three curves agree to better than 5% up to 4 times above threshold. 
Furthermore, all other cases covered in this study fall within the same range. This 
result shows that the influenced of various focusing configurations is mostly accounted 
for by the dependence of pump threshold Pthr (figures 3,6 and 9) and maximal depletion 

Dmax (figures 1,4 and 7). With these two parameters, the conversion D at any pump 

power level P is predicted to within 5% by scaling to the plane wave formula: 



Here D pw (P') is the pump depletion for plane wave beams at A'=0 as a function 
of the ratio P'=P/Pthr of pump power to threshold power. According to ref. [8], this is 
defined implicitly by giving D pw and P' as functions of the signal amplitude A's (see 
equations 4.a and 4.b): 



D(P/P t hr) = D max D pw (P/Pthr) 



(14) 



Dpw(A's) = sin(A' S L')2, P'= (A S 'L')2/sin(A S 'L')2, L'=L/zO. 



(15) 



The maximal pump depletion for plane wave beams at A'=0 is D 



'pw max 



=1.0. 



3. Intensity and phase dependent effects of focusing 

The dependence of SRO conversion on focusing results from both intensity and 
phase sensitive mechanisms. 

Intensity dependence is given by the fact that tighter focusing leads to higher field 
amplitudes, and thereby to enhanced nonlinear conversion. This effect is 
counterbalanced by the reduction in effective gain length which results from stronger 
beam spreading. The latter effect dominates at very short confocal parameters. These 
combined effects are most clearly seen in figure 7 to 9, for the case zp'=zs'. 

If zp'=zs', the nonlinear polarization at the idler wave length has a radial phase 

dependence which, as a function of propagation distance, is that of a Gaussian with the 
same constant confocal parameter and location of focus as the pump and signal, as long 
as the pump retains its Gaussian form. Therefore, conversion can be phase matched to a 
zero order idler mode. On the other hand, if zp' is different from zs', the radial phase 

dependence of the non-linear polarization does not correspond to that of a zero order 
idler mode, and phase matching is frustrated. Therefore, smaller gain than in the phase 
matched case results, and one expects that higher pump levels are required for threshold 
and depletion. Also, higher transverse idler modes are expected to be generated more 
strongly, leading to a corresponding reduction in maximal depletion D max , since the 

different transverse modes reach maximal power at different levels of signal power. 
These qualitative arguments are supported by the focusing dependence of D max 

and required pump levels discussed earlier. For example, D max in figures 1 and 4 peaks 

at values of zp' close to zs'. Furthermore, for zp'>0.7, tighter focusing of the signal 

beam leads to a decrease in D max , and for zp'>0.9 also to an increase in required pump 

power (figures 2 and 5). This is contrary to what would be expected on the basis of 
intensity dependent effects only, and results form the increasing radial phase mismatch. 
For zp'<0.7, on the other hand, tighter focusing of the signal beam increases both the 

signal amplitude and the degree of radial phase match, and therefore gives higher D max 

and lower required pump power. 

The dependence of higher radial idler mode content on focusing is demonstrated 
in table 1 below, which shows the maximal depletion D max and the relative power in the 

zero order radial mode for kp/kp=0.50, zs -1.0 and zp'=0.4, 0.8 and 1.5. Best phase 



matching occurs for zp'=0.8, which has highest D max and highest content of zero order 

radial modes. For calculation of radial mode amplitudes, the confocal parameter of the 
idler modes was set equal to the pump confocal parameter, and focus located at the 
crystal center. 



pump confocal parameter 1.5 0.8 0.4 

maximal pump depletion 0.92 0.97 0.93 

relative zero-mode strength 0.65 0.70 0.62 

Table 1: Conversion parameters for lq/kp=0.5 and bs/L=1.0. 



4. Summary and conclusion 

In this paper, we presented detailed numerical results for the focusing dependence 
of the maximal pump depletion and of the required pump levels for various degrees of 
pump depletion in SROs. We found that maximal depletion occurs for values of the 
pump confocal parameter slightly less than or equal to the signal confocal parameter. 
The focusing dependence of conversion and required power levels is mostly contained 
in the values of pump threshold Pthr and maximal depletion D max . The relative 

depletion D/D max as a function of relative pump power P/Pthr is equal to the result from 

the plane wave theory to within 5%. In particular, D=0.95D max is achieved at about 2 

times above threshold. 

Both intensity and phase dependent effects are important. If the confocal 
parameters of pump and signal are almost equal, the focusing dependence of maximal 
conversion and threshold is due mostly to the dependence of signal amplitude and 
effective gain length. For larger differences between confocal parameters, the radial 
phase mismatch is the dominant factor, and conversion increases with reduced phase 
mismatch, even if a decrease in signal amplitude is involved. 

These results present a simple quantitative and qualitative account of the main 
effects of diffraction in SROs, which can be used to easily predict and assess the 



performance of focused beam based SROs. The main limitations of our model are: a) 
restriction to radially symmetric beams and zero order pump and signal beams, and b) 
neglect of any nonlinear effects other than parametric conversion and dispersion. In 
particular, thermal effects and resulting radial changes in refractive index may be of 
importance in continuous-wave SROs [13]. 

It is a pleasure to acknowledge interesting discussions about the topics presented 
here with W. Bosenberg-who stimulated our interest in SROs-, R. Byer and M. Fejer. 
Thanks are also due to S. Schiller for inviting us to contribute to this special issue. 
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Figure Captions 



Figure 1 : Maximal pump depletion D max as a function of pump confocal parameter. 

Figure 2: Required pump power for D=D max and D=0.95 D max . 

Figure 3. a: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
b S /L=1.0. 



Figure 3.b: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
bs/L=0.7. 

Figure 3.c: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
bs/L=0.5. 

Figure 3.d: Required relative pump power P/Pthr for D=0.95 D max . 

Figure 4: Maximal pump depletion D max as a function of pump confocal parameter. 

Figure 5: Required pump power for D=D max and D=0.95 D max . 

Figure 6.a: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
bs/L=1.0. 

Figure 6.b: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
b S /L=0.7. 

Figure 6.c: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 at 
b S /L=0.5. 

Figure 6.d: Required relative pump power P/Pthr for D=0.95 D max . 

Figure 7: Maximal pump depletion D max as a function of pump confocal parameter for 
bp=bs. 

Figure 8: Required pump power for D=D max and D=0.95 D max , for bp=bs. 

Figure 9. a: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 for 
bp=bs atki/kp=0.33. 

Figure 9.b: Required pump power for threshold (D=0.0), D=0.2, 0.4, 0.6 and 0.8 for 
b P =bs at lq/kp=0.50. 

Figure 9.c: Required relative pump power P/Pthr for D=0.95 D max . 

Figure 10: Relative pump depletion D/D max as a function of relative pump power 
P/Pthr- 
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Figure 1. Maximal pump depletion 
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Figure 2: Pump power at maximal depletion 
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Figure 3. a: Pump power at constant depletion 
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Figure 3.c: Pump power at constant depletion 
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